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Abstract 


Y 

This  investigation  determines  the  induced  currents  on  a  finite 
sized  array  of  parallel  wires  when  illuminated  by  a  plane  wave  with 
varying  incidence  angles  by  the  method  of  moments.  The  arrays  considered 
are  of  various  spacings,  lengths,  and  widths,  ranging  from  5  to  36  wires 
per  wavelength,  and  one  to  two  wavelengths  in  length.  The  effects  of  a 
ground  plane  parallel  to  the  array,  and  located  one  quarter  wavelength 
away,  is  also  studied. 

The  analysis  is  accomplished  by  the  method  of  moments  using 
piecewise  sinusoidal  expansion  functions  and  Galerkin's  method.  An 
algorithm  is  developed  to  accomplish  the  integration  and  matrix  inversion. 
It  was  written  general  enough  for  the  user  to  specify:  number  of  wires, 
number  of  segments,  wire  length,  wire  diameter,  wire  spacing,  spacing 
above  the  ground  plane,  frequency,  and  magnitude  of  incident  electric 
field.  The  results  illustrate  the  various  effects  that  changing  the  wire 
spacing,  wire  length,  number  of  segments,  and  incidence  angle  have  upon 
the  induced  current.  The  results  are  also  compared  to  the  modified 
physical  optics  approximation.  The  results  of  this  investigation  indicate 
that  the  moment  method  is  accurate  enough  to  produce  very  reasonable 
approximations  of  the  induced  current  for  most  applications. 
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ANALYSIS  OF  THE  INDUCED  CURRENTS  ON  A  SECTION  OF  PARALLEL 


WIRES  IN  FRONT  OF  A  GROUND  PLANE 


I.  Introduction 


1.1  Background 

Considerable  work  has  been  accomplished  by  many  authors  in  the 
analysis  of  the  induced  currents  on  wires.  In  many  problems  of  this 
nature,  the  boundary  condition  that  the  tangential  electric  field  equal 
zero  along  the  wire  surface  must  be  satisfied.  This  approach  will  yield 
a  complex  integral  equation  (Pocklington* s  Integral  Equation).  Classical 
solutions  to  this  integral  equation,  and  therefore  for  the  induced 
currents,  are  tedious  and  solvable  for  only  a  few  simple  wire  geometries. 
Any  wire  geometries  whose  surfaces  are  not  easily  describable  in  a  con¬ 
ventional  coordinate  system  are  generally  unsol vable  by  classical  methods. 

Because  of  the  difficulty  in  obtaining  completely  accurate 
solutions  for  arbitrary  wire  geometries,  justifiable  approximations  are 
often  made  which  yield  good,  approximate  solutions.  One  commonly  used 
approximation  is  a  physical  optics  technique.  Even  the  physical  optics 
technique,  however,  is  limited  to  a  few  basic  wire  geometires.  Most  of 
the  previous  studies  are  limited  to  four  basic  cases:  (1)  a  single  dipole, 
(2)  an  array  of  parallel  dipoles,  (3)  single  wires  of  both  finite  and 
infinite  length,  and  (A)  an  array  of  parallel  wires  of  infinite  length. 
Approximate  solutions  have  been  found  for  the  induced  currents  for  the 
above  cases  when  the  wires  are  assumed  to  be  in  free  space.  A  review  of 
the  articles  relevant  to  this  thesis  is  given  in  the  following  section. 


This  thesis  will  be  an  extension  of  case  (4)  above;  the  problem 
to  be  considered  is  an  array  of  this,  infinitely  conducting,  parallel 
wires  of  finite  length.  Approximate  solutions  will  be  found  for  the 
induced  current  when  the  array  is  completely  in  free  space  and  in  the 
proximity  of  a  ground  plane.  At  present,  no  detailed  analysis  of  this 
problem  exists. 

The  need  to  determine  the  induced  currents  on  an  array  of  parallel 
wires  has  arisen  in  several  situations,  particularly  when  a  wire  array 
antenna  is  being  studied.  A  problem  of  this  nature  of  interest  to  the 
Air  Force  is  a  type  of  Cassegrainian  antenna  with  wire  reflector  surfaces. 
The  general  antenna  geometry  may  be  represented  as  in  Figure  1. 


Figure  1.  Side  view  of  a  Cassegrainian  Antenna  with  wire 
reflector  surfaces. 


In  order  to  determine  the  (far  field)  radiation  pattern,  the 
induced  currents  on  each  of  the  wires  must  be  found.  Because  of  the 
complicated  geometry,  completely  accurate  results  are  realistically 
impossible  to  obtain.  By  making  some  justifiable  approximations,  the 
problem  can  be  reduced  to  a  workable  complexity,  and  approximate  solutions 
for  the  induced  current  can  be  obtained. 

A  fundamental  approximation  to  be  made  concerns  the  geometry  of 
the  wires  and  the  ground  plane.  Working  with  the  curved  ground  plane 
and  wire  surfaces  of  the  antenna  is  obviously  undesirable;  approximating 
these  surfaces  as  flat,  rather  than  curved,  considerably  simplifies  the 
problem.  This  approximation  can  be  justified  by  considering  the  entire 
curved  rear  reflector  to  be  subdivided  into  a  large  number  of  perfectly 
flat  sections.  The  approximation  will  approach  the  true  condition  as  the 
number  of  divisions  increases.  Therefore,  the  accuracy  of  the  approximate 
solution  can  be  increased  as  desired  by  simply  dividing  the  reflector 
into  a  greater  number  of  sections.  If  the  induced  current  can  be  found 
for  any  arbitrary  flat  section,  then  by  superposition,  the  approximate 
currents  on  the  entire  reflector  will  be  known  to  any  desired  accuracy. 

The  problem,  then,  has  been  reduced  to  determining  the  induced 
currents  on  the  wires  of  an  arbitrary,  flat  section.  In  order  to  remain 
general,  the  solution  must  be  valid  for  an  arbitrary  section  which  lies 
on  the  edge  of  the  reflector.  To  account  for  this  possibility,  the 
currents  may  be  forced  to  zero  on  at  least  one  end  of  the  array.  However, 
a  more  general  problem  will  be  considered  where  the  current  is  forced  to 
zero  at  both  ends  of  the  array. 
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While  no  exact  analytical  solution  for  wires  of  finite  length 
exists,  the  currents  are  often  assumed  to  be  represented  by  a  modified 
physical  optics  approximation.  The  usual  physical  optics  approximation 
is 

J  =  2n  x  H1  (1) 

s 

where  n  is  a  unit  normal  vector  to  the  scattering  surface 
is  the  incident  magnetic  field  vector 
Jg  is  the  surface  current  density 

The  physical  optics  approximation  is  actually  an  approximate  solution  to 
the  surface  current  density  Jg  on  a  large,  flat  conducting  plate.  It 
seems  reasonable,  then,  that  an  array  of  closely  spaced,  parallel  thin 
wires  can  be  approximated  by  a  conducting  plane,  as  in  Figure  2.  The 
approximation  will  become  increasingly  better  as  the  radius  and  the 
spacing  of  the  wires  approach  zero. 

In  order  to  account  for  the  fact  that  the  current  may  only  flow 
along  the  wire  axis,  the  actual  physical  optics  approximation  is  modified 
as 

j  =  e-(2n  x  H1)  (2) 

s 

where  e  is  a  unit  vector  parallel  to  the  wire  axis  as  in  Figure  2. 

1.2  Review  of  Literature 

As  stated  in  the  previous  section,  considerable  work  in  obtaining 
the  induced  currents  on  wire  structures  has  been  accomplished  by  various 
authors.  This  section  will  briefly  review  those  works  which  are  most 
relevant  to  this  thesis. 


Conducting  plane  approximation  for 


The  first  article  to  be  reviewed  was  written  by  Jack  H.  Richmond 
(Ref  1).  It  was  one  of  the  first  articles  written  discussing  digital 
computer  solutions  of  scattering  problems.  The  paper  introduced  a  new 
technique  for  solving  for  the  induced  currents  on  perfectly  conducting 
and  dielectric  bodies.  Before  the  paper  was  published  in  1965,  high 
speed  digital  computers  were  relatively  unavailable.  The  accepted  tech¬ 
nique  for  solving  problems  of  this  nature  prior  to  1965  were  variational 
and  quasi-static  methods. 

Richmond  states  that  considerable  success  for  scatterers  of 
various  shapes  has  been  established;  however,  these  techniques  (varia¬ 
tional  and  quasi-static  methods)  are  limited  to  bodies  which  are  small 
in  comparison  to  the  wavelength,  or  at  most  are  on  the  order  of  one 
wavelength  in  maximum  diameter.  Larger  scatterers  are  handled  with  the 
aid  of  physical  optics,  geometrical  optics,  and  the  theory  of  diffraction. 

Richmond's  paper  introduced  a  technique  which  generated  a  system 
of  linear  equations  by  enforcing  the  boundry  conditions  at  many  points 
within  the  scatterer  or  upon  its  surface.  Then,  with  the  aid  of  a  digital 
computer,  this  system  of  equations  was  solved  to  determine  the  current 
distribution  of  the  surface,  or  the  coefficients  in  the  mode  expansion  for 
the  scattered  field.  The  distant  radiated  pattern  was  then  found  using 
the  approximated  currents. 

This  technique  is  impractical  without  a  computer,  because  a  large 
system  of  equations  must  be  solved  to  provide  a  reasonably  accurate 
solution.  Richmond  expects  the  linear  equation  solution  to  be  accurate 
for  bodies  of  arbitrary  material,  size,  and  shape.  His  paper  briefly 
reviews  recent  progress  of  the  technique,  discusses  methods  for  reducing 


the  computational  effort,  and  illustrates  an  example  by  considering  a 
plane  wave  to  be  incident  upon  a  single  straight  wire  or  dielectric  rod 
of  finite  length. 

Richmond  begins  by  stating  that  the  electric  field  intensity  is 
the  sum  of  the  incident  and  scattered  intensities.  Since  it  is  known 
that  the  tangential  electric  field  at  the  wire  surface  must  be  zero 
(the  wire  was  assumed  to  be  infinitely  conducting) ,  the  scattered  field 
equals  the  negative  incident  field.  He  writes  the  scattered  field  as 
integral  of  the  current  density  multiplied  by  an  expansion  function  over 
the  length  of  the  wire.  The  expansion  function  used  here  is  a  Fourier 
Series. 

Richmond  states  that  even  if  the  true  current  distribution  were 
actually  known,  the  integral  describing  the  radiated  electric  field 
could  not  be  evaluated  analytically,  except  in  the  form  of  an  infinite 
series.  He  claims,  though,  that  numerical  integration  is  both  possible 
and  feasible  with  a  digital  computer.  It  should  be  noted  that  any 
solution  employing  this  technique  will  only  be  approximate,  because  only 
a  finite  number  of  terms  may  be  realistically  considered. 

Results  are  given  for  the  Fourier  coefficients  of  the  current 
expansion  on  a  wire  illuminated  by  an  incident  plane  wave.  The  length 
of  the  wire  is  varied  from  0.1  to  0.7  wavelengths;  the  diameter  of  the 
wire  is  0.01  wavelengths.  The  results  show  good  agreement  with  the 
experimental  results  (within  10%)  when  the  integrals  are  evaluated  with 
fifth  order  Newton-Cotes  formulas  (which  are  exact  for  fifth  order 
polynomials)  using  1000  terms. 
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Richmond  also  applied  this  linear  equation  technique  to  a 
dielectric  rod  of  0.0625  inches  in  radius,  0.5  inches  in  length.  Results 
include  a  graph  of  the  rod  length  versus  the  echo  area  per  square  wave¬ 
length.  His  results  show  that  the  field  scattered  from  the  rod  differ 
significantly  from  that  of  the  physical  optics  approximation  (by  about 
10%),  but  agree  well  with  experimental  data  (within  10%). 

Therefore,  it  can  be  concluded  that  (1)  the  physical  optics 
approximation  is  probably  not  applicable  to  single  dielectric  rods  and 
wires  of  small  radius,  and  (2)  that  this  linear  equation  technique  using 
infinite  series  expansions  will  not  be  practical  when  arrays  of  wires 
are  considered,  because  Richmond  used  1000  terms  to  obtain  good  results 
for  a  single  wire. 

The  second  article  to  be  reviewed  was  written  by  K.  K.  Mei; 
numerical  solutions  of  a  dipole  antenna  are  considered  (Ref  2) .  Mei 
begins  by  citing  the  integral  equation  used  by  Pocklington  for  z  directed 
dipoles 

_  2  _ j 

/  J(z)[ -  G(z,z')  +  32G(z,z')]dz'  =  -jtoeE  (3) 

L  9z2  z 

where  J(z)  is  the  current  density 

G(z,z^)  is  the  free  space  Green's  function,  - — 

4ttR 

A  numerical  solution  of  the  integral  equation  is  then  obtained  by  approxi¬ 
mating  the  integration  at  a  finite  sum  of  N  points.  A  matrix  equation  is 
therefore  generated.  In  effect,  Mei  has  simply  "relaxed"  the  boundry 
condition  that  the  tangential  electric  field  must  be  everywhere  zero  to 
only  be  zero  at  N  discrete  points.  This  is  known  as  the  point  matching 


solution  of  the  moment  method.  Although  the  point  matching  technique  is 
only  an  approximation,  good  results  were  obtained. 

In  a  separate  article,  Richmond  considered  scattering  by  an 
arbitrary  array  of  parallel  wires  (Ref  3) .  He  considered  three  wires 
not  in  the  same  plane,  30  parallel  wires  in  a  semicircle,  and  15  parallel 
wires  in  an  I  beam  configuration.  Each  wire  is  infinitely  long.  The 
radius  and  spacing  was  0.03  and  0.2  wavelengths,  respectively. 

A  set  of  N  linear  equations  is  generated  by  representing  the 
current  on  each  wire  as  the  unknown  quantity.  This  is  not  a  point 
matching  technique  since  only  one  equation  is  generated  for  each  wire, 
along  the  center  of  each  wire.  Because  each  wire  is  infinitely  long, 
Richmond  writes  the  current  on  each  wire  as  an  expansion  of  a  Hankel 
function.  Thus,  the  electric  field  is  defined  as  the  product  of  a  con¬ 
stant,  a  current  magnitude,  and  a  Hankel  function  expansion. 

This  technique,  then,  is  to  represent  the  current  on  each  of  the 
infinitely  long  wires  in  the  form  of  Hankel  functions.  A  linear  equation 
is  then  written  for  each  of  the  N  wires.  These  linear  equations  are  then 
solved  using  matrix  algebra. 

Results  were  obtained  for  each  of  the  above  wire  configurations. 

A  graph  of  electric  field  vs.  incidence  angle  is  presented  which  compares 
Richmond's  results  with  the  physical  optics  approximation  for  conducting 
strips  of  similar  geometry.  Richmond's  results  compare  favorably.  Thus, 
it  is  evident  that  solutions  for  the  currents  on  an  array  of  wires  can  be 
obtained  by  assuming  that  the  current  is  uniform  along  each  wire.  How¬ 
ever,  this  necessitates  the  wires  to  be  infinitely  long.  Therefore,  this 
technique  is  not  applicable  to  the  problem  of  finitely  sized  arrays. 


9 


In  the  fourth  article  under  review.  Wait  investigates  a  similar 
problem  to  that  of  Richmond,  i.e.,  reflection  from  a  parallel  wire  grid 
(Ref  4) .  Wait  outlines  a  solution  for  the  problem  of  a  plane  wave  which 
is  obliquely  incident  upon  an  infinitely  long  parallel  wire  array.  He 
assumes  that  the  wires  are  small  compared  with  the  separation  and  the 
wavelength. 

The  plane  wave  has  an  electric  field  which  can  be  given  by 


E(x,y,z)  =  Aexp{jB(xcos^ocos9  +  ysintJ^cosQ  -  zsinG)} 


as  shown  in  Figure  3. 


Figure  3.  Wait's  planar  wire  grid  with  incident  plane  wave, 


With  this  assumption,  the  radius  of  the  wires  must  be  small  com¬ 


pared  to  the  wavelength,  since  axial  symmetry  is  assumed.  Because  the 
wires  are  infinitely  long,  the  solution  can  be  found  by  using  Hankel 
functions,  in  a  similar  manner  to  Richmond.  Although  Wait  gives  no  con¬ 
crete  examples  or  experimental  data,  his  theory  agrees  with  results 
obtained  by  Richmond. 

Scattering  by  an  array  of  infinitely  long  parallel  wires  is  also 
the  subject  of  the  fifth  paper  reviewed.  Michel,  Pauchard,  and  Vidal 
describe  a  mathematical  solution  to  the  scattering  problem  by  an  infi¬ 
nitely  long  array  (Ref  5) .  They  claim  that  their  technique  applies  to 
conductors  which  are  (1)  continuous,  (2)  loaded  with  localized  impedances, 
or  (3)  cut  in  equal  length  segments  and  equally  spaced  in  a  colinear 
allignment  of  dipoles.  The  continuous  case  is  most  similar  to  this  thesis. 

They  assume  the  wires  to  be  parallel  and  infinitely  long. 

Equations  are  generated  using  Hankel  functions  which  are  very  similar  to 
those  found  by  Richmond  and  Wait.  The  current  on  each  of  the  conductors 
is  represented  by  the  product  of  constants,  a  current  magnitude,  and  a 
Hankel  function  expansion.  Their  paper,  however,  is  more  general  than 
Wait's  or  Richmond's,  since  theory  accounting  for  the  possibility  that 
localized  impedances  may  lie  on  the  conductors  is  included. 

As  an  example,  they  give  a  double  planar  array  of  conductors 
which  is  alternately  continuous  with  no  impedances,  and  also  continuous 
with  impedances.  These  results  are  compared  to  the  experimental  radia¬ 
tion  pattern  of  a  double  planar  array  illuminated  by  a  uniform,  equi- 
phase  plane  wave.  Their  theoretical  results  agree  quite  well  with 
experimental  data. 
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Thus,  to  this  point,  the  only  geometries  of  wires  which  have  been 
analyzed  in  detail  are  (1)  single  wires  of  infinite  length,  (2)  dipoles, 
or  (3)  wire  arrays  of  infinite  length.  No  author  has  considered  the 
problem  of  a  finite  sized  array  of  closely  spaced,  parallel  wires. 

Richmond,  however,  in  a  subsequent  article,  did  consider  a  wire  grid 
model,  where  the  wires  are  finite  length  (Ref  6). 

Richmond  shows  that  a  point  matching  solution  can  be  developed 
for  scattering  by  conductors  of  arbitrary  shape.  In  his  paper,  he  describes 
a  practical  technique  for  calculating  scattered  fields  of  wire  loops,  con¬ 
ducting  plates,  and  bodies  of  arbitrary  shape.  He  models  a  conducting 
plate  with  intersecting  wire  arrays  as  shuwn  in  Figure  4. 


Figure  4.  Richmond's  wire  grid  model. 


Richmond's  wire  grid  model  is  not  like  the  problem  under  con¬ 
sideration,  since  it  contains  two  intersecting  wire  arrays.  His  case 
does  assume,  however,  that  the  wires  are  of  finite  length.  The  major 
distinction  is  that  the  current  is  modeled  to  flow  along  either  wire  axis 
or  any  combination  thereof.  This  thesis  is  not  considering  a  finite  con¬ 
ducting  plate  but  rather  a  finite  array  of  wires,  where  the  current  is 
allowed  to  flow  only  along  the  wire  axis. 

Richmond  solves  for  the  induced  currents  on  his  wire  grid  model 
by  a  point  matching  technique.  Each  of  the  finite  length  wires  in  the 
wire  grid  model  is  divided  into  an  equal  number  of  segments  of  lengths. 

Then,  the  tangential  electric  field  is  forced  to  zero  only  at  a  single 
point  at  the  midpoint  of  each  segment.  Therefore,  a  system  of  linear 
equations  is  generated  which  can  be  solved  by  the  aid  of  a  digital  com¬ 
puter. 

In  a  later  paper,  which  is  on  the  same  topic,  Richmond  gives  a 
computer  program  for  solving  thin  wire  scatterers  using  sinusoidal  bases 
and  Galerkin  method  (Ref  7) .  As  examples  to  illustrate  the  validity  of 
this  point  matching  technique,  Richmond  gives  results  comparing  experimental 
data  with  point  matching  solutions  for  circular  wire  loops,  square  wire 
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loops,  square  wire  grid  models,  conducting  hemispheres,  and  conducting 
spheres.  Each  shows  good  agreement  (5  to  10%  difference)  with  experi¬ 
mental  results. 

Tne  work  which  presented  much  of  the  basic  moment  method  theory 
used  in  this  thesis  also  merits  a  brief  review.  Stutzman  and  Thiele  pre¬ 
sented  a  good,  fundamental  explanation  of  the  moment  method  in  their  text. 
Antenna  Theory  and  Design  (Ref  8:306-332).  Although  they  omit  much 


detail,  a  derivation  of  Pocklington' s  Integral  Equation  for  a  single  wire 
case  using  the  Lorentz  gauge  condition  and  free  space  scalar  Green’s 
function  is  presented.  They  include  a  section  explaining  how  Pocklington' s 
Integral  Equation  resembles  Kirchoff's  Network  Equations  when  expansion 
functions  are  used.  Two  examples  with  results  are  included:  they  are 
(1)  point  matching  on  a  short  dipole,  and  (2)  a  pulse-pulse  Galerkin 
solution  for  a  short  dipole.  Although  no  examples  are  included  for  the 
piecewise  sinusoidal  expansions,  a  brief  description  of  their  efficiency 
is  given.  Included  in  the  Appendix  is  a  program  for  calculating  the 
currents  on  a  single  array  of  parallel  dipoles.  Because  their  matrix  is 
block  Toeplitz,  it  cannot  be  applied  to  this  thesis,  i.e.,  an  array  over 
a  ground  plane,  since  this  case  does  not  produce  a  block  Toeplitz  matrix. 

The  final  article  reviewed  was  recently  written  (July  ’82)  by 
Kastner  and  Mittra  (Ref  9:673-679).  The  paper  is  an  analysis  of  a 
corrugated  surface  twist  reflector  by  a  spectral  iteration  technique.  It 
is  included  in  this  review  to  illustrate  that  alternative  techniques  for 
analyzing  twist  reflectors  do  exist.  Although  their  corrugated  surface 
twist  polarizer  is  not  exactly  the  same  as  the  Cassegrainian  Antenna 
reflector  (Figure  1)  of  this  study,  both  antennas  are  functionally  similar, 
i.e.,  they  are  designed  to  rotate  the  incident  wave  90°.  Instead  of  using 
a  wire  array  over  a  ground  plane,  their  corrugated  surface  twist  reflector 
uses  a  set  of  thin  metallic  strips  on  top  of  a  dielectric  substrate  which 
is  A/4  thick  and  is  backed  by  a  ground  plane,  thus  forming  a  series  of 
troughs  which  are  A/4  deep. 

They  claim  their  spectral  iteration  approach  is  more  efficient 
than  the  conventional  moment  method  or  mode-matching  techniques,  since 


the  spectral  iteration  approach  requires  no  matrix  inversion.  Without 
including  all  the  details,  the  spectral  iteration  approach  for  this 
problem  is  based  on  the  fact  that  Green's  function  for  both  the  regions 
interior  and  exterior  to  the  trough  is  expressible  in  terms  of  a  Fourier 
type  series.  Because  of  this,  the  authors  indicate  that  in  each  of  these 
regions  the  integral  operator  that  relates  the  E  and  H  fields  is  exactly 
invertible. 

There  are  two  different  iteration  sequences,  depending  upon 
whether  the  incident  wave  is  TE  or  TM.  The  general  sequence,  however,  is 
to  (1)  evaluate  the  outer  region  operator  via  two-step  fast  Fourier 
transform  (FFT) ,  (2)  obtain  H  or  E,  (3)  evaluate  the  inner  region  oper¬ 
ator  and  estimate  E  of  H,  and  (4)  repeate  as  desired  for  accuracy.  The 
reflection  coefficient  is  computed  at  each  step  of  the  iteration,  and 
convergence  is  attained  with  three  or  four  iterations. 

The  authors  claim  that  the  use  of  the  FFT  procedure,  combined  with 
the  small  number  of  iterations  needed  to  attain  convergence,  result  in 
a  considerable  amount  of  saved  computation  time  as  compared  to  the  moment 
method  and  mode-matching  techniques  which  require  a  large  matrix  inversion 
Results  are  included  for  a  sample  twist  reflector,  and  the  authors  cite 
the  results  to  be  "very  good". 

Thus,  the  existing  situation  is  that  a  general  technique  (the 
moment  method)  exists  for  approximating  the  induced  currents  on  finite 
sized  arrays.  Richmond  has  demonstrated  that  the  moment  method  produces 
accurate  results  for  his  finite  sized  wire  grid  model  and  other  geometries 
The  problem  of  determining  the  induced  currents  on  an  array  of  finite 


sized  wires,  though,  has  not  been  solved  to  date.  The.  most  similar 
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geometries  to  the  problem  of  this  thesis  which  have  previously  been 
solved  are  (1)  an  array  of  parallel  wires  of  infinite  length,  and 

(2)  wires  of  finite  length  arranged  in  two  intersecting,  perpendicular 
arrays  as  in  Richmond's  wire  grid  model,  Figure  4. 

1.3  Problem 

The  problem  under  consideration  is  to  determine  the  induced 
currents  on  an  array  of  parallel  wires  of  finite  length,  both  in  free 
space  and  in  front  of  a  ground  plane.  This  thesis  will  compare  these 
results  with  those  of  the  modified  physical  optics  approximation. 

1.4  Scope/Assumptions 

This  thesis  will  determine  the  induced  currents  on  an  arbitrary 
section  (array  of  parallel  wires  in  front  of  a  ground  plane) .  The  wires 
will  be  of  finite  length,  and  will  have  a  small  but  finite  diameter. 
Length,  diameter,  and  spacing  will  be  varied  in  order  to  gain  some 
insight  of  their  effect  on  the  currents.  As  initial  values,  the  diameter 
will  be  approximately  0.1  of  the  wire  spacing;  the  spacing  will  be  25 
wires  per  wavelength. 

To  make  the  problem  solvable  under  given  constants,  several 
assumptions  must  be  made.  This  thesis  will  assume  the  following: 

(1)  Ground  plane  is  perfectly  flat  and  perfectly 
conducting  (o  =  °°) . 

(2)  Wires  are  perfectly  conducting  (a  =  °°) . 

(3)  Wires  are  0.25  wavelength  in  front  of  the 
ground  plane. 


(4)  The  incident  plane  wave  is  at  an  arbitrary 

angle  of  incident  with  matched  polarization. 


(5)  The  wires  have  a  small  enough  diameter  and 
are  far  enough  apart  so  that  there  is  no 
variation  of  current  around  the  circumference, 
only  along  the  length. 


1 . 5  Approach 

The  analysis  will  be  broken  down  into  two  cases.  The  first  case 
will  be  the  analysis  of  the  currents  on  a  section  of  wire  with  no  ground 
plane  present.  The  second  case  will  be  the  analysis  of  a  section  of  wires 
with  a  ground  plane  present.  Image  theory  will  be  used  to  model  the  wires 
and  ground  plane  as  two  rows  of  parallel  wires  separated  by  0.5  wave¬ 
lengths.  In  both  cases,  the  polarization  will  be  assumed  to  be  parallel 

to  the  wire  axis  (TM  mode) . 

z 

The  actual  analysis  of  the  currents  in  each  case  will  be  accom¬ 
plished  by  using  the  method  of  moments,  as  described  by  Stutzman  and 
Thiele.  The  method  of  moments  is  a  technique  for  solving  an  integral 
equation  which  can  be  readily  implemented  on  a  computer.  This  will  be 
done  by  approximating  an  integral  equation  by  a  set  of  linear  algebraic 
equations  with  the  induced  current  being  the  unknown. 

Each  wire  will  be  divided  into  N  segments,  with  each  segment 
having  an  unknown  current  function.  This  current  function  can  be  any 
type  of  continuous  functions,  i.e.,  pulse,  triangle,  etc.  If  pulse 
functions  are  used,  the  technique  is  known  as  point  matching,  since  only 
the  center  point  of  each  segment  is  forced  to  obey  the  equation.  For 
this  study,  sinusoidal  functions  will  be  assumed.  By  also  assuming  that 
both  the  source  segments  and  testing  functions  are  sinusoidal,  the  pro¬ 
cedure  is  referred  to  as  Galerlcin  Method.  The  reason  for  using  sinusoidal 
bases  (rather  than  other  functions)  is  that  the  current  is  assumed  to  be 
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sinusoidal  in  nature,  and  therefore  more  accurate  solutions  should  be 
obtainable  with  fewer  of  segments,  thus  reducing  the  size  of  the  matrix 
to  be  inverted.  In  this  way,  the  numerical  computational  effort  by  the 
computer  is  reduced. 

Since  the  incident  field  is  assumed  to  be  known,  we  have  a  matrix 

equation  with  one  unknown  of  the  form  V  =  Z  I  ,  wiiere 

V.  is  the  known  Incident  Field, 

a 

Z^  is  an  impedence  matrix  representing  the  array  of  wires, 
each  of  which  is  divided  into  N  segments. 

I  is  the  unknown  current  distribution. 
a 

I  can  then  be  found  by  using  the  computer  to  perform  a  matrix  inversion 
of  Z.  Finally,  these  results  will  be  compared  to  the  modified  physical 


optics  approximation. 


II.  Development  of  Theory 


<?• 


2.1  Wire  Geometry 

As  stated  in  Chapter  I,  this  thesis  will  be  a  study  of  an  array 
of  parallel  wires,  both  in  free  space  and  in  the  proximity  of  a  ground 
plane.  In  both  cases,  the  relative  distance  between  any  two  arbitrary 
segments  must  be  known,  because  the  integral  equation  is  a  function  of 
this  distance.  In  order  to  determine  this  distance,  it  is  necessary  to 
define  an  orientation  of  the  array  on  a  coordinate  system.  Although  any 
coordinate  system  may  be  used,  a  rectangular  coordinate  system  is  chosen 
here  for  convenience. 

2.1.1  Free  Space  Geometry 

For  this  case,  the  array  of  wires  is  considered  planar,  and  can 
be  expressed  in  two  dimensions  (temporarily  ignoring  the  wire  radius) . 

The  array  will  be  defined  as  in  Figure  5,  where 
N  =  number  of  segments  per  wire 
M  =  number  of  wires 
WSY  =  spacing  between  adjacent  wires 
WL  =  wire  length 
L  =  segment  length 

In  this  manner,  any  arbitrary  segment  can  be  specified  by  an  array  coordi¬ 
nate  (n,m).  Because  there  are  M  wires,  each  with  N  segments,  there  will 
be  a  total  of 

Q  =  MN  (5) 

segments . 


19 


The  ground  plane  is  assumed  to  be  perfectly  conducting.  When  the 


array  is  located  near  the  ground  plane,  the  total  electric  field  at  a 
point  P  will  be  the  sum  of  two  components:  (1)  the  direct  field,  and 
(2)  the  reflected  field.  Using  Image  Theory,  the  array  and  ground  plane 
can  be  modelled  as  two  parallel  arrays  separated  by  twice  the  original 
separation,  as  in  Figure  6. 


GROUND  I MAGS 

ARRAY  PLANE  ARRAY 


<r-~  0.2^  -4- 0.25 A— 


Figure  6.  Image  Theory  model  for  an  array 
near  a  ground  plane. 

Therefore,  since  the  array  is  0.25  wavelengths  above  the  ground  plane, 
the  equivalent  geometry  is  two  parallel  arrays  separated  by  0.5  wave¬ 
lengths.  The  orientation  for  this  case  will  be  similar  to  that  in 
Figure  5,  except  that  the  image  array  will  exist  in  the  X  =  -0.5A  plane. 


For  this  case,  the  segments  are  relabelled  so  that  the  total  number  of 
segments 

Q  =  2MN  (6] 

for  easier  implementation  to  the  computer  program  (see  Appendix  A) . 


2.1.3  Distance  Between  Two  Arbitrary  Segments 

The  purpose  of  this  section  is  to  define  the  geometry  for  calcu¬ 
lating  the  distance  between  two  segments.  For  purpose  of  identification, 
one  segment  is  labelled  "observation"  segment  (obs) ,  with  its  location 
with  reference  to  the  origin  given  by 

robs  =  +  yy  +  zz  (7) 

Similarly,  the  other  segment  is  labelled  "source",  with  its  location  with 
reference  to  the  origin  given  by 


r 

source 


>A  ,  ^  A  ,  >  A 

=  xx  +  yy  +  zz 


(8) 


For  notational  purposes,  primed  variables  will  refer  to  source  segments, 
and  unprimed  variables  will  refer  to  observation  segments. 

The  distance  between  obs  and  source  segments,  then,  is  the  magni¬ 
tude  of  the  difference  between  these  two  vectors,  hence 


R  =  I  r 


-  r  , 

source  obs' 


/(x"-x)2  +  (y'-y)2  +  (z"-z)2 


(9) 


as  in  Figure  7. 


22 


Figure  7.  Representation  of  separation  between  segments. 


2.2  Fields  At  An  Arbitrary  Wire  Segment 


This  section  develops  the  theory  for  the  electric  field  at  an 
arbitrary  observation  segment  caused  by  a  current  on  an  arbitrary  source 
segment.  Beginning  with  the  time-harmonic  form  of  Maxweel's  Equations 


V  X  E  =  -jtoyH 


V  X  H  =  jwe'E  +  3 


and  using  the  definition  for  the  magnetic  vector  potential  A, 


H  =  V  X  A 


Eq  (11)  can  be  written  as: 


V  X  H  =  V  X  V  X  A  =  jre'E  +  j 


Using  the  vector  identity 


V  X  V  X  A  =  V(V-A)  -  V2A  (14) 

Eq  (13)  is  also  equal  to 

V(V,A)  -  V2A  4  ju)C>(-jwyA  -  V4>)  +  J  (15) 

because  the  scalar  potential  $  is  defined  from: 

E  4  -jwyA  -  V$  (16) 

Rewriting  Eq  (15)  as 

V2A  +  oa2ye^A  -  V(jcoe'<i>  +  V*A)  =  -j  (17) 

and  employing  the  Lorentz  condition  (Ref  8: 10)  to  specify  the  divergence 
of  the  vector  potential 

V*A  =  (18) 

then  Eq  (17)  reduces  to: 

V2 A  +  w2ye'A  =  -J  (19) 


which  is  known  as  the  "vector  wave  equation". 

A  solution  to  this  vector  wave  equation  is  (Ref  8:13) 

A-  "1 

source  4ttR 
volume 


(20) 


or  in  terms  of  the  scalar  Green's  function 


1 


ft  O' 


where  iJj(R)  = 


,-jkR 

4ttR 


A  =  ^  Jip  (R)  dv^ 

source 

volume 


,  Green's  function. 


(21) 


From  the  definition  of  the  scalar  potential,  Eq  (16),  E  can  be  expressed 
in  terms  of  A.  This  can  be  done  by  taking  the  divergence  of  Eq  (18) 


V(V* A)  =  -jwe'V$ 


(22) 


and  solving  for  V$,  yielding 


V$  =  -- 


j  we 


■  V(V*A) 


(23) 


Substituting  Eq  (23)  into  Eq  (16)  leaves 


E  =  -jWUoA  +  V (7*A) 


juie 


(A +  VZA) 


(24) 


or,  because  the  propagation  constant  B2  =  w2pe'  ,  then 


E  =  - 


(BZA  +  V2A) 


jwe 


(25) 


By  finally  inserting  the  Green's  function  wave  equation  solution  Eq  (20) 
into  Eq  (25),  the  desired  result  is  obtained: 


E  =  [V2  ip  ( R)  +  B2^(R)]Jdv> 

jwe  source  r  y 

volume 


(26) 
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Some  clarification  is  needed,  however,  on  what  is  meant  by 
"source  volume".  If  a  typical  source  segment  is  defined  in  cylindrical 
coordinates  as  in  Figure  8, 


then: 


Figure  8.  A  typical  segment  in  free  space. 


.  L  2tt  a 

E  =  Iff  .  72iKR)  +  B2^(R)]jp"dp"d(J)"dz' 

^  o  o  o 


(27) 


It  will  he  assumed  here  that  the  conductivity  of  air  approximately  equals 
the  conductivity  of  free  space,  so  that  z'  =  z  .  Because  it  was  assumed 
that  the  conductivity  of  the  wire  was  infinite  (a  =  °°) ,  all  of  the  current 
density  is  assumed  to  exist  only  upon  the  wire  surface  (at  p  =  a) ,  not  on 
the  interior  region.  Hence, 


L  2tt 


jwe 


/  /  [V?<KR)  +  B2tKR>Jdd<Tdz' 


(28) 


o  o  o 
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Since  it  was  stipulated  that  the  wire  must  be  very  small  in 
diameter,  it  can  be  assumed  that  the  current  only  flows  along  the  wire 
axis,  not  transverse  to  it.  Therefore,  this  assumption  implies  that  only 

A  A 

a  Z  electric  field  will  exist;  the  Z  component  of  the  electric  field 
expressed  in  Eq  (28)  is  then: 

_  L  2ir  -2 

Ez  =  tot  1  f  W  +  ^2^(r)  IJdfy'dz'  (29) 

J  o  o  o  z 

Also  because  the  wires  are  of  small  diameter,  it  can  be  assumed  that  the 
current  density  is  equally  distributed  on  the  circumference  of  the  wire 
surface  in  an  infinitely  thin  shell  (only  at  p  =  a) ,  so 

I(z')  =  J(z")2TTa  (30) 

Then,  substituting  Eq  (30)  into  Eq  (29)  and  integrating  yields  the  final 
desired  result: 

Ez  =  /  [yr  >KR)  +  B2iKR)]l(z/)dz'  (31) 

J  o  o  z 


4 


I 


Before  continuing,  it  should  be  noted  that  the  assumption  of  the 
current  existing  only  upon  the  wire  surface  will  necessitate  a  modifica¬ 
tion  to  Eq  (9).  The  observation  points  of  each  segment  are  along  the 
wire  axis;  with  the  above  assumption,  the  source  points  are  on  the  wire 
surface.  This  difference,  the  wire  radius  a,  must  be  taken  into  consid¬ 
eration.  If  the  radius  were  not  included,  it  would  be  possible  for  R  to 
equal  zero  when  the  source  and  observation  segments  were  the  same.  If  R 
did  equal  zero,  Green's  function  would  become  infinite,  which  would 
require  consideration  of  the  shape  of  the  principle  volume  and  the 
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principle  value  of  the  integral.  Fortunately,  this  is  not  necessary  in 
this  case.  Eq  (9)  should  therefore  be  modified  to  include  the  radius  as: 


R  =  /a2  +  (x'-x)2  +  (y^-y)2  +  (z'-z)2 


2.3  Method  of  Moments 

The  method  of  moments  is  a  technique  for  numerically  approxi¬ 
mating  a  solution  t^  an  integral  equation  of  the  first  kind,  like  Eq  (31). 
Schelkunoff  demonstrated  that  this  technique  is  analogous  to  solving 
Kirchoff's  network  equation  (Ref  8): 


Z  Z  ,1  ,  =  V 

-  j  qq  q  q 


where 


Q  =  total  number  of  segments 


q  =  q  source  segment 


q'  =  q  observation  segment 
By  defining  a  function  K(R)  as 


K(R)  =  ~ -  [■“  <K R)  +  62i^(R)] 


then  Eq  (31)  can  be  expressed  as: 


E  -  /  K(R)I(Odz' 
z 

o 


Furthermore,  the  current  on  a  wire  can  be  approximated  by  a  series  of 
orthogonal  expansion  functions  such  that 


KO  ~  E  I  F  (z') 
n=l  n  n 
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will  be  a  system  of  Q  equations 


Ez(zp  K  IjfCzj.zp  +  I^Cz^zp  +  ...  +  Ipf  (z  ^ ,  Zq> 

Ez(z2)  ~  I1f(z2>zp  +  I2f(z2,z2)  +  ...  +  Igf (z2 J zq) 


Ez(zq)  *  IjfCzp.z')  +  I2f(zQ,z2)  +  ...  +  Iqf(zQ,z') 


(41) 


which  can  be  more  conveniently  expressed  in  matrix  notation  as: 


[V  ]  =  [Z  .]  [I  J 

q  qq  q 


(42) 


The  I^'s  can  then  be  found  by  simply  premultiplying  Eq  (42)  by  [z  ..] 
obtaining : 


[i  A  =  [Z  A  1  [V  ]  (43) 

q  qq  q 

The  previous  few  paragraphs  are  an  explanation  of  the  general 
form  of  the  method  of  moments  (Ref  8:310-312).  Specific  forms  of  the 
method  of  moments  are  obtained  by  choosing  different  expansion  functions 
in  Eq  (36)  . 

Although  any  set  of  orthogonal  functions  may  be  used  for  expan¬ 
sion  functions,  the  simplest  and  most  intuitively  obvious  choice  may  be  a 
pulse  function: 

1  for  z'  in  L 

F  (z')  = 

n  0  else  (44) 


where  L  is  a  segment  of  the  wire  length  WL/N. 


This  choice  of  expansion  functions  is  commonly  used  with  the  "point 
matching  solution".  It  is  so  labelled  because  the  boundry  condition 
(tangential  electric  field  equals  zero)  is  only  enforced  at  N  discrete 
points.  A  physical  interpretation  of  this  is  shown  in  Figure  9  (Ref  8:311). 


Figure  9.  Staircase  approximation  to  an  actual 
current  distribution. 

Exactly  accurate  solutions  are  obtained  using  this  method  if  N 
equals  infinity.  In  practice,  however,  good  approximations  are  obtained 
by  making  N  sufficiently  large.  N  cannot  be  made  arbitrarily  large, 
however,  due  to  difficulties  encountered  with  the  large  computational 
effort . 

Stutzman  and  Thiele  state  that  approximately  fifty  segments  per 
wavelength  are  needed  to  obtain  reasonable  results  (Ref  8:324).  For  a 
single  wire  as  they  studied,  the  point  matching  approach  is  adequate. 
However,  when  a  large  array  is  encountered  (as  in  this  thesis  of  25  wires 


per  wavelength) ,  it  becomes  evident  that  the  point  matching  approach 
will  become  computationally  unfeasible.  For  example,  if  this  thesis  were 
to  attempt  a  point  matching  approach,  the  size  of  the  array  to  be  inverted 
would  be  (50  segments  times  25  wires  implies)  a  1250  by  1250  array;  this 
is  clearly  an  unreasonable  task  for  all  but  the  most  advanced  computer 
systems.  Obviously,  a  different  expansion  function  must  be  chosen. 

This  thesis  will  use  piecewise  sinusoidal  expansion  functions 
(PWSEF)  defined  as 


which  are  pictured  in  Figure  10  (Ref  8:324). 


Figure  10.  Piecewise  sinusoidal  expansion  function. 


Each  segment  will  overlap  with  adjacent  segments  to  approximate  the 
current  on  the  entire  wire,  as  shown  in  Figure  11. 


Figure  11.  Five  segment  piecewise  sinusoidal 
current  approximation. 

PWSEF  are  chosen  because  the  actual  current  is  anticipated  to  be 
sinusoidal,  and  therefore  the  closest  approximation  should  be  PWSEF. 

The  greatest  advantage  in  using  PWSEF  instead  of  pulse  functions  is  the 
inherent  reduction  in  segments.  To  obtain  the  same  degree  of  accuracy, 
Stutzman  and  Thiele  find  that  "almost  ten  times  fewer  segments  are 
required  for  the  PWSEF  as  for  the  pulse  expansion  function"  (Ref  8:325). 

Therefore,  using  PWSEF,  only  about  five  segments  per  wavelength 
are  needed.  This  reduction,  then,  makes  the  study  of  a  25  wire  per  wave¬ 
length  array  feasible  using  the  moment  method  with  PWSEF.  Now,  only  a 
125  by  125  array  (250  by  250  for  ground  plane  case  because  of  image  array) 
need  be  inverted.  This  is  possible  for  many  computer  systems. 


In  order  to  obtain  the  best  possible  approximation,  the  boundary 


condition  (tangential  electric  field  equal  zero)  must  be  enforced  as 
"closely"  as  possible.  The  point  matching  technique  is  not  an  exception¬ 
ally  good  approximation  because  the  boundary  condition  is  only  enforced 
at  one  point  in  each  segment.  Instead,  a  better  technique  is  to  attempt 
to  enforce  the  boundary  condition  in  an  "average"  sense  (as  many  points 
as  possible)  over  the  segment.  This  is  what  is  done  in  an  approach  known 
as  the  method  of  weighted  residuals  (Ref  10). 

T6  S 

A  residual  field  E^.  can  be  defined  as 

tan 


pres 

Etan 


+  E 


tan 


(46) 


Clearly,  the  exact  solution  occurs  when  E  =  0  .  In  any  approximate 

1TGS 

solution,  Etan  ^.0  ;  however,  the  best  approximation  would  be  to  make 
E  approach  zero  in  some  "average"  sense.  This  can  be  done  by  defining 

t  all 

a  Weighting  Function  W^(z)  such  that: 

L 

/  W  (z)E^eS(z)dz  =  0  n= 1 , 2 , . . . N  (47) 

n  tan 
o 


The  weighting  function  W^(z)  is  similar  to  the  expansion  function 
F^(z)  ;  it  can  be  any  set  of  orthogonal  functions.  If  the  weighting  and 
expansion  functions  are  chosen  to  be  the  same  (in  this  case,  piece-wise 
sinusoidal),  then  the  technique  is  known  as  the  Galerkin  method.  As 
previously  stated,  this  thesis  is  using  the  Galerkin  method;  therefore, 
W^(z)  =  F^(z)  •  Then  using  Eq  (46),  Eq  (47)  is  equivalent  to: 

L  I. 

/  W  (z)E®  (z)dz  +  f  W  (z)E*  (z)dz  =  0  (48) 


Since  a  current  on  a  wire  will  produce  a  scattered  field  E  ,  and 
because  there  are  only  z  directed  wires,  Eq  (31)  can  be  interpreted  to 
represent  the  scattered  field  which  was  caused  by  the  induced  current. 

Eq  (31)  can  then  be  substituted  into  Eq  (48)  for  E^an  ,  yielding: 


/  Wn(z)  [v  !'  /  <KR)  +  32tKR)]l(z')dz']dz 
L  o  o 


-/  W  (z)E1  dz 
L  n  tan 


The  right  side  of  Eq  (49)  is  recognized  to  be  the  "average" 
induced  voltage,  ;  it  will  be  discussed  in  the  following  section.  The 
left  side  is  equal  to  the  induced  current  multiplied  by  the  cross  imped¬ 
ance  Z  ..  between  the  source  and  observation  segments.  It  can  be 
qq 

simplified  by  realizing  that: 


■h  *(R)  -  -  ■si’ 


JL  «R>  =  jl-.j  *(R) 


(50a) 


(50b) 


Using  Eas  (50a)  and  (50b)  in  Eq  (31)  and  integrating  twice  by 
parts  yields  (Ref  8:329): 


Ez  -  »<R>  b  I(*'> + 1(2')  £ 

o  z 


+  Hz')  +  62i(z')]'MR)dz'  (51) 

^  O  o  ^ 
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Since  the  current  is  being  approximated  by  PWSEF,  for  a  single 
til 

source  segment,  Eq  (51)  becomes: 


E  =  — 
z  toe 


[*<*>  +  Fg'(z')  ♦<*>! 


+  &  V(2'>  +  BJrq.(2')]<KR)dz'  (52) 

O  Li 


An  interesting  simplification  occurs  because  PWSEF  were  chosen.  Using 
Eq  (44)  and  twice  differentiating,  it  is  discovered  that: 


7^2  F  '(*')  = 


-32F  a*') 


This  indicates  that  the  bracketed  expression  in  the  integrand  of  Eq  (52) 
is  zero  (Ref  9:370).  Eq  (52)  then  reduces  to: 


sh  [«R>  sr  V<*'>  +  V<z'>  £  «'<R)] 

o  n  source 

segment 


Evaluating  Eq  (54)  over  a  typical  source  segment  (pictured  in 
Figure  10),  and  denoting  the  distance  from  the  obs  point  to  » 

and  Zq'+i  as  R  ,  R  -  »  and  R  >  respectively,  Eq  (54)  equals: 

E  =  -J-  [Ui(R  )  F  ,(z  ,)  +  F  Az  ,)  i|)(R  ) 

z  U)Cq  q  9z  q  q  q  q  Sz^  q' 


-  ST  W-l1  •  W-P  i  ^V-P 


+  ♦‘V+P  ST  VV+P  +  V<V+1>  A  ^V+P 


<|i(R  .)  -  F  .(;:  0  -  F  ,(::  ,)  .J  \' (E  )] 

q  dz  q  q  q  q  Jz  q 
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Then,  the  expression  for  Z  *  can  be  given  using  Eq  (58): 


O' 


Z  -  =  /  I  E  dz 


qq 


obs 

seq 


q  Z 


-  C  / 


2 

z  sin8(z  -  z  ,)  q+1  sin8(z  -  z) 

,.q _ q~  1  t _ q+1 _ 


sin8(z  -  z  . ) 
z  .  q  q-1 

q-1 


+  / 


sinB(z  -  z  )' 
zq  q+1  q' 


[  j  30  (- 


-jBR  -  , 

2  J  q-1 


R  >  ,sin8(z  -  z  -  , ) 
q-1  q  q-17 


cos8(z  s , ,  -  z  ^ )  -jBR  . 
2  q+1  q  e  JM  q 


sin6(zq.+i 


z  >)  R  - 

q  q 


e  j^Vjfl 


R  ^,.sin8(z  -  z  ) 

q+1  q+1  q' 


]dz 


(59) 


Stutzman  and  Thiele  present  a  similar  but  less  detailed  derivation 
for  their  simpler  case  of  a  single  wire.  Eq  (59)  agrees  with  their 
result  if  appropriately  simplified  (Ref  8:331). 


2.4  Induced  Voltage  by  Incident  Plane  Wave 

This  section  describes  the  q1"  voltage  matrix  element,  Vq  .  From 
Eq  (49),  the  "average"  weighted  voltage  is: 


V 

q 


=  /  W  (z) 


E  dz 
tan 


(60) 


Etan  must  be  determined  for  a  plane  wave  of  arbitrary  incidence  angle 
with  matched  polarization,  as  stated  in  the  initial  assumptions.  An 
arbitrary  plane  wave  can  be  specified  by  the  "plane  wave  solution"  to  the 
vector  wave  equation  as: 


E  = 


E 

o 


+ j  Bn • r 


(61) 
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E  is  a  constant  vector,  to  be  specified 
o 

8  is  the  propagation  constant 

n  is  the  propagation  vector 

r  is  the  radius  vector  from  the  origin 

Since  matched  polarization  is  assumed,  then  =  Eqz  ,  because 
the  wires  are  z  directed.  By  expressing  the  radius  vector  r  in  terms  of 
rectangular  coordinates,  then 


E  =  E  z  e+j£((sinecos<£)x  + 
o 


(sin0sin<f>)y  4-  (cos0)z) 


so  that: 


(6: 


Zq  Sin8(z  -  z  ,)  Zq+1  sinB (z  ,  -  z) 

v  _  _T  r  _ _  q  f  j.  r  _ _ 1 

q  sinB (z  -  z  ,)  sinB(z  .  -  z  )J 

zq_l  <1  q-1  zq  q+1  *J 


E  z  e+jS((sin9cos4>)x  +  (sin0sin<j))y  +  (cos0)z))dz 
o 


(6 


Now,  both  Z  ^  and  V  have  been  expressed  in  a  form  suitable  for 

qq  q 


computer  implementation  (see  Appendix  A),  Results  for  different  wire 
spacings,  lengths,  and  incidence  angles  are  given  in  the  following 


chapter. 


III.  Results 


This  chapter  presents  the  results  of  the  computer  program  and 
compares  them  with  the  physical  optics  approximation.  The  parameters 
investigated  include:  (1)  wire  spacings,  (2)  number  of  segments  per 
wire,  (3)  wire  lengths,  and  (4)  varying  incidence  angles.  Unless  other¬ 
wise  indicated,  all  results  were  obtained  with  the  same  set  of  input 
values.  These  are: 

Number  of  wires  =  as  indicated 
Segments  per  wire  =  as  indicated 
Wire  length  =  1.0  meters 
Wire  radius  =  0.005  meters 

X  wire  spacing  =  0.50  meters 

Y  wire  spacing  =  determined  by  — mTr-Wr— - 

'  a  J  width  of  array 

Frequency  =  300. OMH^ 

Propagation  constant  6  =  2tt/A 

Ie  I  =  1.0  volts 
1  o 1 

Width  of  array  -  1.0  A  =  1.0  meters 

The  current  is  in  units  of  AMPS,  voltage  in  VOLTS,  and  length  in  METERS. 

For  each  of  the  different  cases  studied,  both  current  density 
magnitude  and  phase  is  presented  in  either  a  "horizontal  cut"  or  a 
"vertical  cut".  These  terms  need  to  be  properly  defined  for  accurate 
interpretation.  A  horizontal  cut  is  the  current  on  the  center  segment  of 
each  wire  across  the  width  of  the  array.  For  example,  with  5  segments 
per  wire,  a  horizontal  cut  presents  the  3rd,  8th,  13th, . . . ( Q—  2 ) th  seg¬ 
ments.  Conversely,  a  vertical  cut  presents  the  current  of  the  center  of 


the  wire  array  for  the  length  of  the  wire.  For  example,  with  N=5  and  M=7, 
a  vertical  cut  presents  the  current  on  the  center  wire,  segments  16,  17, 
18,  19,  and  20. 

3. 1  Wire  Spacings 

A  fundamental  question  to  be  considered  is  for  what  range  of 
wire  spacings  is  this  moment  method  technique  applicable.  It  is  known 
that  the  technique  can  be  used  successfully  for  a  single  wire  case,  as 
demonstrated  by  Stutzman  and  Thiele  (Ref  8).  However,  an  objective  of 
this  thesis  was  to  determine  if  a  parallel  wire  array  could  sufficiently 
model  a  flat,  conducting  plane.  This  involves  using  "enough"  wires  to 
achieve  a  "reasonable"  approximation. 

As  was  stated  in  Chapter  I,  it  is  expected  that  the  moment  method 
results  should  become  more  accurate  as  the  number  of  wires/wavelength  is 
increased.  Results  in  Figures  12  and  13  verify  this  expectation  to  a 
certain  extent.  The  purpose  of  establishing  these  results  is  to  deter¬ 
mine  if  a  trend  exists  as  the  number  of  wires  is  increased. 

For  a  normally  incident  plane  wave  upon  a  flat,  conducting  sheet, 
the  physical  optics  approximation  would  yield  a  current  density  of; 


i1  =  Ze+jkx 


=  Zi 


+jkx 


376.7 


J  =  e*2nxH 
s 


376.7 


5. 31mA/m 


This  does  not  totally  agree,  however,  with  what  the  actual  distribution 
is  expected  to  be.  One  would  expect  the  actual  current  to  be  basically 


uniform  (flat)  across  the  center  of  the  sheet,  with  perhaps  a  perturbation 


Figure  13.  Current  on  a  vertical  cut  of  a  1A2  array  with  various 
wire  spacings  for  a  normally  incident  plane  wave. 


near  the  edge  caused  by  the  abrupt  boundary  change.  The  current  along 
each  wire  (vertical  cut)  should  be  sinusoidal,  with  the  endpoints  of 
each  wire  equal  to  zero. 

It  is  evident  from  Figures  12  and  13  that  the  moment  method  seems 
to  be  more  accurate  than  the  physical  optics  approximation.  The  hori¬ 
zontal  cut  indicates  a  basically  flat  distribution  across  the  center  with 
a  perturbation  at  the  edge,  as  is  physically  expected.  The  physical 
optics  approximation  is  flat  and  of  the  correct  order  of  magnitude,  but 
does  not  represent  any  effects  caused  by  the  edges.  The  vertical  cut 
is  symmetric  and  sinusoidal,  as  expected.  The  current  is  zero  at  the 
ends  of  the  wire,  but  Figure  14  does  not  apparently  indicate  this 
because  the  current  at  the  center  of  each  sinusoidal  piece  was  plotted, 
not  the  current  at  the  exact  endpoint  of  the  wire.  Note  that  the  physical 
optics  approximation  does  not  show  the  currents  approaching  zero  at  the 
ends  of  the  wire. 

The  current  distributions  also  tend  to  "flatten  out"  in  the 
center  region  as  the  number  of  wires  is  increased,  as  expected.  A  word 
of  caution  is  in  order,  however.  Although  this  flattening  trend  is 
present  through  36  wires/wavelength  (which  should  be  enough  for  most 
applications) ,  it  may  not  continue  to  exist  for  arbitrarily  large  numbers 
of  wircs/wavelength .  It  was  initially  assumed  that  the  wire  diameter 
was  "small  enough"  so  that  there  was  no  circumferential  current  variation 
on  the  wires.  If  the  spacing/diameter  ratio  becomes  too  small,  this 
assumption  may  become  invalid.  Although  the  lower  limit  of  this  ratio 
has  not  been  identified  by  this  thesis,  good  results  are  achievable  as 
low  as  ( 1/36) /0. 005  =  5.55:1.  This  point  may  be  an  aspect  of  future  study 
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As  presented  in  the  theory,  more  accurate  results  should  be 
obtained  by  using  a  greater  number  of  segments/wire.  Theoretically, 
there  is  no  upper  limit;  there  is  an  upper  limit  imposed,  however,  by 
the  increased  computational  effort  incurred  and  the  available  computer 
resources.  The  question  this  section  investigates  is  whether  increasing 
the  number  of  segments  (beyond  five)  is  worth  the  additional  effort. 

Figures  14  and  15  display  the  induced  current  (horizontal  and 
vertical  cuts,  respectively)  for  a  IX2  array  with  five  wires.  By 
increasing  the  number  of  segments,  one  would  expect  a  more  accurate  dis¬ 
tribution  of  the  current  along  the  length  of  the  wire,  with  little  or  no 
change  across  the  array  (in  the  horizontal  cut) .  Note  that  the  physical 
optics  approximation  is  unchanged  for  each  case,  since  changing  the  number 
of  segments  has  no  meaning  with  this  approximation. 

Results  show  good  agreement  with  actual  expected  behavior.  The 
horizontal  current  distribution  changes  little  in  both  general  shape 
and  magnitude  for  all  five  cases.  This  is  reasonable,  since  one  would 
not  expect  the  current  across  the  wires  to  be  appreciably  affected  by 
changing  the  sampling  increment  along  the  length  of  the  wire.  The  verti¬ 
cal  current  also  changes  as  expected.  The  midpoint  current  remains 
relatively  constant  while  the  endpoints  decrease  nearly  in  order  of 
magnitude . 

This  sixfold  increase  in  the  number  of  segments  shows  no  signifi¬ 
cant  change  in  the  horizontal  current  or  the  center  region  of  the  vertical 
current,  but  docs  indicate  the  expected  decrease  near  the  endpoints  of  the 
wire.  Since  it  is  known  that  the  current  must  be  zero  at  the  ends  of  each 
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Figure  14.  Current  on  a  horizontal  cut  of  a  1AZ  array  with  various 
wire  segment  divisions  for  a  normally  incident  plane  wave. 


Figure  15.  Current  on  a  vertical  cut  of  a  1A2  array  with  various 
wire  segment  divisions  for  a  normally  incident  plane  wave. 
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wire,  the  region  of  interest  is  the  relatively  uniform  center  region. 
Because  increasing  the  number  of  segments/wavelength  beyond  five  con¬ 
tributes  little  additional  information,  it  can  be  concluded  that  five 
segments/wavelength  is  sufficient  to  obtain  reasonably  accurate  results. 
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3.3  Wire  Length 

Figures  16  and  17  present  the  horizontal  and  vertical  current 
distributions  on  a  IX  wide  array  for  four  values  of  wire  lengths.  The 
purpose  of  this  section  is  to  determine  if  the  moment  method  will  reflect 
the  expected  changes  along  the  wire  length. 

Results  for  the  horizontal  cut  are  reasonable.  All  four  curves 
are  symmetrical  with  respect  to  the  center  wire,  and  tend  to  flatten  out 
as  the  wire  length  is  increased.  This  is  reasonable,  since  as  the  array 
size  is  increased.,  one  would  expect  the  edges  to  have  less  of  a  distorting 
effect  upon  the  center  region,  thus  allowing  the  current  to  become  more 
uniform.  The  decreasing  trend  is  explained  by  considering  the  vertical 
cut.  Figure  17. 

Currents  on  the  vertical  cut  also  change  appropriately  as  the 
wire  length  is  increased.  Figure  17  indicates  a  trend  to  form  two  main 
lobes  (d)  as  the  length  is  increased  to  2. OX,  which  seems  quite  acceptable. 
As  these  two  lobes  are  formed,  the  center  segment  (it 10)  value  must  be 
decreased  accordingly;  this  in  turn  causes  the  decrease  in  the  horizontal 
cut.  Figure  16.  Therefore,  the  primary  effect  of  increasing  the  wire 
length  is  to  cause  the  appropriate  number  of  lobes  to  form  along  the 
length  of  the  wire,  and  to  cause  the  horizontal  current  to  shift  accord¬ 
ingly.  Note  that  the  physical  optics  approximation  reflects  no  change 
in  either  cut  as  the  wire  length  is  increased. 
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Figure  16.  Current  on  a  horizontal  cut  of  a  IX  wide  array  with  various 
wire  lengths  for  a  normally  wide  plane  wave. 


Current  on  a  vertical  cut  of  a  IX  wide  array  with  various 
for  a  normally  incident  plane  wave. 


3.4  Varying  Incidence  Angle  Effects 


The  horizontal  and  vertical  current  distributions  are  determined 
for  different  sized  arrays  when  the  plane  wave  incidence  angles  (0,  <j>) 
are  varied  along  their  major  axis,  i.e.,  when  0  is  varied  and  4>=0,  and 
when  <J>  is  varied  and  when  8=tt/2.  Phase  plots  of  each  are  also  included 
in  this  section,  since  there  will  be  phase  variations  in  both  cuts. 

The  first  array  considered  is  IX2,  using  26  wires  of  five  segments 
each.  Results  for  this  example  are  presented  in  Figures  18  through  23. 

For  0  variations.  Figures  18  and  20  present  the  horizontal  and  vertical 
currents,  respectively.  Physically,  one  would  expect  little  change  in 
shape  for  the  horizontal  cut  as  0  is  varied,  since  all  segments  will  still 
lie  on  the  same  phase  front  for  all  6's.  This  is  exactly  what  is  shown 
by  Figure  18.  For  this  and  all  other  examples,  the  physical  optics 
approximation  of  the  current  and  phase  is  calculated  by: 


with 

then 


inQ) £  £ejB((sin0cos4>)x  +  (sinOsincfOy  +  (cos0)z) 

2EQ(sin0)  (cost})) 

_____ 
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and  ?  J  =  -y-  (xsinBcost})  +  ysin0sin<}>  +  zcos0) 
S  A 


(64) 
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For  each  0,  the  physical  optics  approximation  is  of  correct  order  of 
magnitude,  but  does  not  represent  the  distortion  of  the  current  near  the 
edges . 

The  currents  for  the  $  variations,  Figures  19  and  21,  are  also 
in  agreement  with  expected  results.  The  horizontal  current  is  symmetric 
at  normal  incidence,  and  varies  as  <p  is  varied.  Note  that  the  physical 
optics  approximation  shows  only  a  shift  in  magnitude,  which  does  not 
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gure  18.  Current  on  a  horizontal  cut  of  a  26  wire 
segment  array  for  0  varying,  <J)  =  0. 


gure  19.  Current  on  a  horizontal  cut  of  a  26  wire 
segment  array  for  <f»  varying,  9  =  tt/2. 


ure  20.  Current  on  a  vertical  cut  of  a  26  wire, 
egment  array  for  0  varying,  cf>  =  0. 


Current  on  a  vertical  cut  of  a  26  wire 


CURRENT  PHASE  (iBGREES)  ' 


Figure  22.  Horizontal  cut  of  current  phase  on  a 
26  wire,  5  segment  array  for  various  incidence  angles 


Vertical  cut  of  current  phase  on  a 
segment  array  for  various  incidence  angles. 


seem  to  agree  with  actual  physical  behavior.  The  vertical  cut  of  the 
current.  Figure  21,  however,  is  symmetric  for  each  case.  This  is 
entirely  reasonable,  since  a  change  in  the  <J>  component  of  the  incident 
plane  wave  would  not  be  expected  to  affect  the  shape  of  the  current  in 
the  0  direction.  There  is  only  a  shift  in  magnitude,  which  is  probably 
caused  by  the  change  in  the  horizontal  current. 

The  phase  relationships  for  both  cuts  are  shown  in  Figure  23  and 
24.  As  expected,  the  phase  in  the  horizontal  cut  is  dependent  upon  <p, 
but  independent  of  0.  The  physical  optics  approximation  predicts  a 
linear  phase  change  for  each  cf>  as  shown,  which  is  calculated  from  Eq  (64) . 
The  moment  method  technique  also  indicates  a  similar  linear  change,  with 
a  slight  distortion  near  the  edges  as  shown.  The  vertical  cut  of  the 
phase.  Figure  23,  is  also  reasonable.  As  expected,  there  is  no  change 
in  phase  for  $  variations,  only  0  variations.  The  physical  optics 
approximation  predicts  a  linear  change,  as  shown.  The  moment  method 
results  show  a  similar,  but  not  exactly  linear  variation,  which  may  be 
caused  by  reflections  from  the  ends  of  the  wire,  which  will  tend  to 
induce  a  "standing  wave"  on  the  wire  as  0  is  decreased. 

A  second  case  was  also  studied  to  investigate  current  behavior 
on  an  array  for  various  incidence  angles.  Its  dimensions  were  chosen  as 
IX  wide  by  2A  long  (13  wires  of  10  segments  each)  to  determine  the 
currents,  especially  in  the  vertical  cut.  Since  the  previous  example  has 
already  established  that  the  horizontal  currents  are  independent  of  0, 
and  the  vertical  currents  are  independent  of  4>,  these  graphs  are  omitted. 

The  horizontal  cut,  Figure  24,  indicates  a  definite  change  in 
the  distribution  as  <f>  is  changed,  as  exported.  Again  note  that  the 


Figure  24.  Current  on  a  horizontal  cut  of  a  13  wire 
10  segment  (2A  long)  array  for  $  varying,  9  =  7t/2. 


physical  optics  approximation  reflects  no  change  in  shape,  only  a  shift 
in  magnitude.  Also  of  interest  is  the  fact  that  the  relative  changes 
in  the  horizontal  cut  for  the  moment  method  results  are  less  than  the 
corresponding  results  for  the  previous  example  (see  Figure  19) .  Again, 
this  is  an  indication  that  the  center  region  becomes  more  "isolated" 
from  current  variations  as  the  size  of  the  array  is  increased. 

The  vertical  cut  for  9  variations,  Figure  25,  is  also  reasonable. 
The  current  for  normal  incidence  illumination  is  symmetric,  and  has  the 
two-lobed  structure,  as  expected  for  a  2\  long  wire.  As  0  is  varied, 
the  current  becomes  asymmetric,  and  in  fact  resembles  a  standing  wave 
pattern.  The  physical  optics  approximation  shows  no  change  in  shape  for 
different  0's,  only  a  shift  in  magnitude. 

The  phase  relationships  for  this  example  are  shown  in  Figures  26 
and  27.  The  horizontal  cut  of  the  phase  agrees  well  with  physically 
expected  behavior.  For  changes  in  9,  there  is  no  change  in  general  shape, 
only  a  shift  in  magnitude,  as  expected.  As  0  is  decreased,  the  rate  of 
phase  change  across  the  array  tends  to  increase,  indicating  that  the  rate 
of  phase  change  is  proportional  to  the  incidence  angle.  The  vertical  cut 
of  the  phase.  Figure  27,  shows  similar  agreement.  Changes  in  4*  have 
little  effect  on  the  vertical  phase,  except  a  small  shift  in  magnitude, 
as  expected.  As  the  incidence  angle  increasingly  deviates  from  normal, 
the  rate  of  phase  change  increases  proportionally.  Note  that  the  physical 
optics  approximations  predicts  an  exactly  linear  phase  shift,  whereas  the 
moment  method  predicts  the  same  general  trend,  but  also  includes  minor 
variations  and  distortions  near  the  edges/ends. 


Vertical  cut  of  current  phase  on  a  13  wire, 
.2).  long)  array  for  various  incidence  angles 


3.5  Array  0.25 A  Above  a  Ground  Plane 


The  current  magnitudes  and  phases  for  an  18  wire,  5  segment 
array  0.25A  above  a  ground  plane  are  presented  in  Figures  28  through  31. 
Again,  because  of  the  independence  of  6  and  horizontal  currents  and  <}>  and 
vertical  currents,  these  graphs  are  omitted. 

The  horizontal  cut  with  varying  <j>  is  shown  in  Figure  28.  As 
expected,  the  normal  incidence  case  is  flat  and  symmetric,  with  increasing 
changes  as  <j>  is  increased.  The  interesting  point  to  consider  is  how 
these  results  differ  from  a  same  sized  array  not  over  a  ground  plane,  as 
Figure  19.  The  physical  optics  approximation  for  the  two  situations 
(both  over  and  not  over  a  ground  plane)  are  the  same.  This  is  because 
the  array  is  being  modelled  as  a  completely  reflecting  ground  plane  itself 
which  in  effect  completely  shields  the  actual  ground  plane  behind  it  from 
the  incident  wave  (at  least  for  normal  incidence — this  is  not  true  for 
all  incidence  angles,  and  this  point  will  be  elaborated  upon  subsequently) 
The  moment  method  technique  shows  that  the  current  in  both  situations 
have  the  same  general  magnitude,  but  that  the  current  on  the  horizontal 
cut  of  the  array  above  a  ground  plane  increases  at  a  greater  rate  across 
the  array.  This  is  reasonable,  as  one  would  expect  the  ground  plane  to 
reflect  the  incident  wave  back  onto  the  array  and  increase  the  current 
density. 

The  vertical  current  cut,  Figure  29,  is  also  somewhat  different 
than  its  free  space  counterpart,  Figure  20.  Again,  the  physical  optics 
approximation  is  the  same  for  both  situations.  The  moment  method 
results,  on  the  other  hand,  are  of  the  same  general  order  of  magnitude, 
but  are  different  in  shape  for  each  situation.  Here  again,  the  moment 


Figure  28.  Horizontal  cut  of  current  density  on  an  18  wire 
5  segment  array  above  a  ground  plane  for  <f>  varying,  0  =  ■ 
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gure  29.  Vertical  cut  of  current  density  on  an  18  wire 
segment  array  hX  above  a  ground  plane  for  9  varying,  <fi  ! 


ure  30.  Horizontal  cut  of  current  phase  on  an  18  wire, 
egment  array  above  a  ground  plane  for  various  incidence 


method  results  appear  more  accurate  than  the  physical  optics  approxima¬ 
tion,  since  different  results  would  be  expected  for  different  situations. 

Similar  conclusions  can  be  made  concerning  the  horizontal  and 
vertical  phase  relationships  for  both  situations.  The  horizontal  and 
vertical  phases  for  the  array  over  a  ground  plane  (Figures  30  and  31) 
are  expected  to  differ  from  the  array  in  free  space  results  (Figures  22 
and  23).  Again,  there  is  no  difference  between  each  for  the  physical 
optics  approximation.  The  moment  method  technique  differ,  as  expected. 

An  important  point  of  the  geometry  of  the  problem  should  be 
noted.  As  was  previously  stated,  the  physical  optics  approximation 
models  the  array  as  a  conducting  sheet,  and  therefore,  the  actual  ground 
plane  is  completely  shielded  from  the  incident  wave  at  normal  incidence. 
At  any  oblique  angle,  however,  part  of  the  incident  wave  will  pass  by 
the  edge  of  the  array,  and  still  illuminate  the  ground  plane;  this  in 
turn  may  then  be  reflected  back  onto  a  different  section  of  the  array 
(see  Figure  32).  This  effect  is  not  accounted  for  in  the  results  of  the 
physical  optics  approximation  as  presented  in  the  previous  example.  To 
properly  account  for  this  effect  would  drastically  complicate  the 
presently  simple  approximation,  and  would  therefore  defeat  the  whole 
purpose  of  using  it.  The  moment  method,  on  the  other  hand,  has  already 
taken  this  effect  into  account  by  virtue  of  Image  Theory,  with  no  addi¬ 
tional  modification  or  complexity. 


gure  32.  Top  view  of  geometry  illustrating 
flection  effects. 


IV.  Conclusions  and  Recommendations 


4.1  Conclusions 

The  objective  of  this  thesis  was  to  determine  the  induced  currents 
on  an  array  of  wires  both  in  free  space  and  in  front  of  a  ground  plane 
using  the  method  of  moments,  and  then  to  compare  the  results  to  those 
obtained  by  the  modified  physical  optics  approximation.  The  induced  cur¬ 
rents  for  the  moment  method  were  determined  by  developing  a  general  algo¬ 
rithm  (Appendix  A) .  Results  were  obtained  for  a  variety  of  different  cases. 

In  general,  moment  method  results  agree  well  with  expected  physi¬ 
cal  behavior,  but  do  not  agree  closely  with  the  physical  optics  approxi¬ 
mation  for  most  cases.  The  moment  method  and  physical  optics  approxima¬ 
tions  tend  to  agree  in  general  magnitudes,  but  sometimes  differ  signifi¬ 
cantly  with  respect  to  the  general  shape  of  the  current  distribution.  For 
example,  consider  an  array  illuminated  by  a  normally  incident  plane  wave. 

The  physical  optics  approximation  would  indicate  the  current  to  be  exactly 
uniform  across  the  entire  array;  the  moment  method  shows  relative  unifor¬ 
mity  in  the  center  region,  with  a  perturbation  near  the  edges.  For  this 
case,  the  moment  method  appears  to  be  more  accurate,  since  one  would 
physically  expect  a  deviation  near  the  edge  caused  by  the  abrupt  boundary 
change. 

Also,  the  physical  optics  approximation  shows  no  change  in  the 
shape  of  the  distribution  as  the  incidence  angles  are  varied;  it  only 
reflects  a  change  in  magnitude.  This  docs  not  appear  to  be  physically 
realistic.  It  can  be  concluded  that  the  physical  optics  approximation 
does  net  adequately  represent  "edge"  effects.  The  physical  optics 


approximation  also  does  not  reflect  the  fact  that  the  current  must 
equal  zero  at  the  endpoints  of  each  wire.  It  does,  however,  produce  a 
reasonably  close  approximation  to  the  current  along  the  center  region  of 
the  wires.  From  this,  it  can  be  generally  concluded  that  the  physical 
optics  approximation  may  be  useful  for  larger  sized  arrays,  to  obtain  a 
"rough"  idea  of  the  current  magnitude,  but  not  useful  to  determine  the 
shape  of  the  distribution  with  any  great  degree  of  accuracy.  Its 
greatest  advantage  is  in  its  simplicity.  For  a  large  enough  array,  it 
provides  a  good  idea  of  the  current  magnitude  with  minimal  effort. 

The  moment  method,  on  the  other  hand,  requires  considerably  more 
effort  (i.e.,  large  matrix  inversions),  but  appears  to  be  far  more 
accurate.  The  accuracy  can  be  increased  by  increasing  the  number  of 
wires  and  segments,  but  this  cannot  be  done  arbitrarily,  since  there  will 
$  be  an  upper  bound  imposed  by  the  available  computer  resources.  However, 

reasonably  accurate  approximations  of  the  currents  are  obtainable  with 
a  modest  number  of  segments,  as  is  evidenced  by  this  thesis.  Some  other 
conclusions  of  the  moment  method  technique  are: 

1.  The  accuracy  of  the  solution  tends  to  increase  as  the  wire 
spacing  is  decreased,  i.e.,  as  more  wires/wavelength  are  used.  Reasonable 
results  are  obtainable  for  as  few  as  five  wires/wavelength,  but  as  many 
as  30  wi  1  provide  an  obviously  better  result.  The  reader  is  cautioned 
that  the  wire  spacing/diameter  ratio  may  not  be  decreased  too  small, 
since  doing  so  may  invalidate  the  assumption  that  there  is  no  circum¬ 
ferential  current  variation  on  the  wires.  Although  this  thesis  made  no 
definite  determination  of  the  lower  limit  of  this  ratio,  good  results 
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were  obtained  with  a  ratio  as  small  as  5.55:1.  The  actual  lower  limit 
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may  in  fact  be  substantially  less  than  this  value. 

2.  The  accuracy  of  the  solution  will  increase  as  the  number  of 
segments/wire  increases.  Although  there  is  theoretically  no  upper  limit, 
little  additional  accuracy  is  gained  when  the  number  of  segments  used 

is  increased  beyond  10/wavelength.  However,  very  reasonable  results  are 
obtainable  with  as  few  as  five  segments/wavelength. 

3.  This  moment  method  technique  appears  to  adequately  represent 
the  various  effects  caused  by  changing  the  wire  lengths  and/or  incidence 
angles  of  the  illuminating  plane  wave.  Results  indicate  that  the  current 
distribution  develops  a  second  lobe  as  expected  when  the  wire  length  is 
increased  to  2. OX.  Also  indicated  was  the  independence,  of  the  horizontal 
current  to  changes  in  0,  as  well  as  the  vertical  current  to  changes  in  <j>. 
The  moment  method  also  indicates  the  distortion  near  the  ends  of  the 
array,  the  trend  toward  zero  current  at  the  endpoints  of  the  wire,  and 
the  different  distributions  of  current  for  the  array  in  free  space  as 
opposed  to  the  array  above  a  ground  plane.  The  physical  optics  approxi¬ 
mation  reflects  none  of  these. 

Thus,  the  moment  method  works  well  in  approximating  the  currents 
on  arrays  both  in  free  space  and  near  a  ground  plane.  The  physical 
optics  approximation  provides  a  very  "rough"  idea  of  the  current  magni¬ 
tude  at  the.  center  of  the  array,  but  does  not  reflect  many  of  the 
important  aspects  of  the  shape  of  the  distribution,  especially  near  the. 
edges/ endpoints  of  the  array.  The  physical  optics  approximation  may  be 
more  useful  for  larger  sized  arrays,  however,  where  the  center  region 
is  more  independent  of  "edge  effects”.  For  smaller  arrays  (on  the  order 
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of  the  ones  in  this  thesis),  however,  the  physical  optics  approximation 
does  not  provide  a  very  accurate  solution. 

4.2  Recommendations 

Although  very  good  results  can  be  obtained  with  the  values  used 
in  this  thesis,  even  greater  accuracy  may  be  desired.  To  accomplish 
this,  one  may  choose  to  incorporate  a  larger  number  of  wires/wavelength. 
As  previously  stated,  this  may  invalidate  the  assumption  that  there  is 
no  circumferential  current  variation  on  the  wires.  Therefore,  a  logical 
extension  of  this  work  could  be  to  better  determine  the  lower  limit  of 
the  wire  spacing/diameter  ratio.  Or,  perhaps  even  better  if  possible,  to 
exclude  this  assumption  entirely  and  rework  the  theory  to  account  for  a 
possible  circumferential  current  variation. 

Secondly,,  the  major  disadvantage  of  the  moment  method  is  the 
large,  time-consuming  matrix  inversion  required.  A  paper  was  recently 
published  (July  1982)  by  Kastner  and  Mittra  (Ref  9:673-679)  describing 
a  new  spectral  iteration  technique  which  analyzes  a  similar  problem  with 
no  large  matrix  inversions;  instead,  it  employs  a  Fourier-series 
approach.  Possible  future  study  may  include  a  detailed  comparison  of  the 
two  techniques. 
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Appendix  A:  Computer  Program 


This  appendix  contains  the  computer  program.  Included  below  is 
a  brief  description  of  the  general  development  and  operation  of  this 
program. 

DEVELOPMENT 

Before  beginning,  it  should  be  noted  that  a  program  (for  a  single 
wire  case)  has  already  been  written  (Ref  8:568).  This  program,  however, 
is  significantly  different;  it  allows  for  the  more  general  case  of  wires 
located  near  a  ground  plane.  The  major  difference  between  the  two  pro¬ 
grams  is  that  this  program  does  not  inherently  assume  that  a  (block) 
Toeplitz  matrix  will  be  obtained.  It  does  not  exploit  the  fact  that  a 
Toeplitz  matrix  is  highly  symmetric  (and  therefore  more  easily  inverted); 
in  fact,  it  assumes  the  most  general  case  that  the  matrix  is  totally 
asymmetric . 

The  program  actually  solves  Eq  (43)  using  Eqs  (59)  and  (63).  In 
order  to  accomplish  the  matrix  inversion,  it  calls  upon  the  International 
Mathematics  and  Statistics  Library  (IMSL)  functions  DCADRE  and  LEQ2C. 

The  IMSL  is  a  set  of  subroutines  which  should  be  available  at  most  com¬ 
puter  installations.  Some  minor  problems  developed  with  the  use  of 
DCADRE,  however,  which  are  worth  mentioning. 

DCADRE  is  the  IMSL  subroutine  to  evaluate  a  real  integral,  given 
the  integrand,  the  limits  (finite),  and  the  maximum  allowable  error, 

AERR.  The  problem  occurred  in  trying  to  use  DCADRE  to  evaluate  the 
complex  integral  in  E.[  (59)  and  (63).  The  solution  used  was  to  break 
the  complex  integral  into  real  and  imaginary  parts,  evaluate  each 


separately,  and  then  recombine  into  complex  form.  Also,  because  of  the 
repetitious  nature  of  the  integrand,  several  common  terms  were  extracted 
and  made  into  function  subprograms,  as  shown  below. 


Zq  sinS(Z  -  Z  ,) 

INTC  =  r  :7W~T  f  - [BRCKTl]dz 

sxn(8*DZ)  ^  sin(8*DZ) 

q-i 


INTD  = 


q+1  sinB(Z  ,  -  Z) 

30  /  ^  . /I  -!— -  [BRCKTl]dz 


sin(8*DZ)  z  sin(B*DZ) 

q 


,n  Zq  sinB(Z  -  Z  ,) 

30  '  [BRCKT2]dz 


INTE  sin(8*DZ)  J  sin(  *DZ) 

q-1 


S  CM 


such  that 


INTF  =  — 


q+1  sin8(Z  -  Z  ) 

30  '  -3-^ - 3_  [BRCKT2]dz 


sin(B*DZ) 


f 


sin(8*DZ) 


Total  Integral  =  (INTC  +  IND)  +  j (INTE  +  INTF) 


where 


BRCKT1  =  {- 


sin(8*R  ^_j) 


2cos(8*DZ) sin(8*R  -) 


R 


q-1 


sin(8*R 

+  — ■■  -a  ■tL) 
I\  >  .  « 

q  +1 
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(A- la) 


(A- lb) 


(A-lc) 


(A- Id) 


(A-2) 


(A- 3a) 


BRCKT2  =  { 


cos(6*R 


R  .  , 
q  -1 


2cos(3*DZ)cos(8*R  ') 


cos(8*R  >.,) 
+  — 

IV  >  ,  « 

q  +1 


(A- 3b) 


and  R  ..  .  ,  R  ,  and  R  ^ ,  are  evaluated  by  functions  RBOT,  RMID, 
q  -1  *  q  q  +1  J 

and  RTOP,  respectively. 


OPERATION 


The  program  was  written  in  FORTRAN  V.  It  requires  the  use  of 
IMSL  subroutines  DCADRE  and  LEQ2C.  Because  of  the  matrix  manipulations, 
it  requires  a  larger  than  average  amount  of  memory  and  CPU  time. 

Although  the  amount  will  vary  depending  upon  the  system  used,  the  user 
should  request  approximately  256K  words  memory  and  9000  seconds  (150  min.) 
of  CPU  time  for  a  250  element  array. 

The  user  should  specify  N,  M,  WSY,  etc  ,  according  to  the  input 
format  statement.  The  output  is  dependent  upon  the  four  control  param¬ 
eters  SWEEP,  PV,  PZ,  and  TYPE. 

SWEEP  must  be  an  integer  value  0  or  1.  0  indicates  that  only  an 

incident  plane  wave  of  normal  incidence  (THETA  =  tt/2,  PHI  =  0)  will  be 
evaluated.  1  indicates  that  the  incidence  angle  will  be  varied  in  ir/16 
increments  along  the  THETA  and  PHI  axis  from  0  to  r/2. 

PV  must  be  an  integer  value  0  or  1.  0  indicates  that  the  voltage 

array  V^(q)  is  not  to  be  printed.  1  indicates  print  V^(q)  . 

PZ  must  be  an  integer  value  0,  1,  or  2.  0  indicates  do  not 

print  impedance  array  Z  ..  .  1  indicates  to  print  only  the  first  row  of 
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Z  .  .  2  indicates  to  print  entire  Z  arrav.  It  is  recommended  that 

qq  qq 

option  0  or  1  be  used,  since  2  can  produce  a  tremendous  amount  of  unneces¬ 
sary  output. 

TYPE  must  be  an  integer  value  0  or  1 .  0  indicates  that  the  free 

space  case  is  to  be  studied,  i.e.,  no  ground  plane.  Here,  Q=MN  so  the 
DIMENSION  statement  is  MN  in  magnitude.  1  indicates  the  ground  plane 
case.  Here,  Q=2MN,  so  each  array  in  the  DIMENSION  statement  should  be 
accordingly  doubled. 

NOTE:  WA  and  WK  arrays  in  the  DIMENSION  statement  are  work  areas 
for  the  IMSL  subroutines.  See  any  IMSL  Operator's  Guide  for  further 
details. 


DATA  OK,  PRINT 
IT  OUT  FOR 
VERIFICATION 


Figure  A- 


Computer  program  flowchart. 
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Figure  A-l  (Continued) 
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Figure  A-2.  Computer  program. 
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